Success-driven distribution of public goods promotes cooperation but preserves defection 
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Established already in the Biblical times, the Matthew effect stands for the fact that in societies rich tend to get 
richer and the potent even more powerful. Here we investigate a game theoretical model describing the evolution 
of cooperation on structured populations where the distribution of public goods is driven by the reproductive 
success of individuals. Phase diagrams reveal that cooperation is promoted irrespective of the uncertainty by 
strategy adoptions and the type of interaction graph, yet the complete dominance of cooperators is elusive due to 
the spontaneous emergence of super-persistent defectors that owe their survival to extremely rare microscopic 
patterns. This indicates that success-driven mechanisms are crucial for effectively harvesting benefits from 
collective actions, but that they may also account for the observed persistence of maladaptive behavior. 

PACS numbers: 89.75.Fb, 87.23.Ge, 89.75.Hc 



The Gospel of St. Matthew states: "For to all those who 
have, more will be given." Roughly two millennia latter 
sociologist Robert K. Merton [ li] was inspired by this writ- 
ing and coined the term "Matthew effect" for explaining dis- 
crepancies in recognition received by eminent scientists and 
unknown researchers for similar work. A few years earlier 
physicist and information scientist Derek J. de Solla Price [2] 
actually observed the same phenomenon when studying the 
network of citations between scientific papers, only that he 
used the phrase cumulative advantage for the description. To 
most physicists, however, preferential attachment will be best 
known due to the seminal paper by Barabasi and Albert f^, 
who used the concept ingeniously to explain the emergence of 
scaling in growing random networks. Other common varia- 
tions of the phrase include rich-get-richer and success-breeds- 
success [4], all implying that initial advantages are often self- 
amplifying and tend to snowball over time. No matter the 
wording, be it the accumulation of wealth [5] or citations fl], 
the making of new friends [7], or the longevity of one's ca- 
reer |8|], this simple yet fascinatingly powerful phenomenon 
arguably influences many facets of our existence. 

In this paper we investigate how the Matthew effect affects 
the evolution of cooperation in the public goods game. The 
public goods game BSfr llOll is played in groups and captures the 
essential social dilemma in that collective and individual in- 
terests are in dissonance. Players must decide simultaneously 
whether they wish to contribute to the common pool or not. 
All the contributions are then multiplied to take into account 
synergetic effects of cooperation, and the resulting amount is 
divided equally among all group members irrespective of their 
initial decision. From the perspective of each individual, de- 
fection is clearly the rational decision to make as it yields the 
highest personal income if compared to other members of the 
group. However, if nobody decides to invest the group fails 
to harvest the benefits of a collective investment and the soci- 
ety evolves towards the "tragedy of the commons" 111 111 



whom to reciprocate with. Social enforcement, on the other 
hand, may work well, although it is challenged by the fact 
that it is costly (see B14I1 for a review). Other prominent ways 
of promoting cooperation in public goods games include the 
introduction of volunteering B15I1 . social diversity by means 
of complex interaction networks JT^], heterogeneous wealth 
distributions 11711 . and institutionalized punishment ll 811 . 

Inspired by the seminal works on games on coevolution- 
ary and social network s 119142311 . the most recent advances on 
this topic II24T - I28I1 (see 112911 for a review), as well as the seem- 
ing omnipresence of the Matthew effect in social interactions, 
we consider the public goods game where the distribution of 
multiplied contributions is not equally shared amongst all the 
group members, but rather it depends on the evolutionary suc- 
cess of each individual. Naturally, the more successful an in- 
dividual is the higher its share of the public good. 

Assuming structured interactions, L 2 players are arranged 
into overlapping groups of size G such that every player is 
surrounded by its G — 1 nearest neighbors. Accordingly, each 
individual belongs to g — G different groups. Initially each 
player on site x is designated either as a cooperator (s x = C) 
or defector (s x = D) with equal probability. Cooperators 
contribute a fixed amount (here considered being equal to 1 
without loss of generality) to the public good while defectors 
contribute nothing. The sum of all contributions in each group 
is multiplied by the factor r > 1 and the resulting public goods 
are distributed amongst all the group members. If s x = C the 
payoff of player x from every group g is Pq = M x rN^ /G—l 
and if s x = D the payoff is P 9 D = M x rN 9 c /G, where 7V£ i s 
the number of cooperators in group g while M x is the fac- 
tor by means of which the Matthew effect is introduced. Ini- 
tially all players have M x — 1, and so without further mod- 
ifications the setup returns the classical spatial public goods 
game l3(ill . Here, however, each time player x successfully 
enforces its strategy on another player y, M x = M x + A and 
The My — My — A, where A > is a free parameter. Thus, in the 



sustenance of cooperation in sizable groups of unrelated indi- 
viduals, as is the case by the public goods game, is particularly 
challenging since group interactions tend to blur the trails of 
those who defect. Unlike by pairwise interactions, reciprocity 
I l_2j, LL3[] often fails as it is not straightforward to determine 



next round player x will receive a higher share of the public 
goods while player y will receive, to the same extent, a smaller 
one. Note that since the fact that player x was able to enforce 
its strategy on player y already implies that the former is more 
successful, this simple coevolutionary rule will strengthen this 
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FIG. 1 : (Color online) Critical R at the D — > C + D phase transition 
in dependence on A for three different values of K. Success-driven 
distribution of public goods strongly decreases the multiplication fac- 
tor needed for cooperation to survive, irrespective of K. Inset depicts 
the probability of strategy adoption w in dependence on the payoff 
difference for the three K [same color (gray scale) as in the main 
panel]. Square lattice with G = 5 was used as the interaction graph. 



further whilst at the same time additionally degrading player 
y, thus concisely introducing the Matthew effect into the pub- 
lic goods game. Importantly, the coevolutionary rule is strat- 
egy independent and at no point in time assumes any public 
goods being lent or not spent, i.e., L~ 2 ^2 x M x = 1 at all 
times. 

Monte Carlo simulations are carried out comprising the 
following elementary steps. A randomly selected player x 
plays the public goods game with its G partners as a mem- 
ber of all the g — 1, . . . , G groups, whereby its overall pay- 
off is thus P Sx = J2 g Pi x ■ Next, player x chooses one of 
its nearest neighbors at random, and the chosen co-player y 
also acquires its payoff P s in the same way. Finally, player 
x enforces its strategy s x onto player y with a probability 
w{s x -> a v ) = 1/{1 + exp[{P Sy - P S J/GK]}, where K 
quantifies the uncertainty by strategy adoptions due to errors 
in decision making or incomplete information (see inset of 
Fig. Q]), and G normalizes the effect for different interaction 
graphs 1 30]. Each Monte Carlo step (MCS) gives a chance for 
every player to enforce its strategy onto one of the neighbors 
once on average. The average frequencies of cooperators ipc) 
and defectors (pd) were determined in the stationary state af- 
ter sufficiently long relaxation times. Depending on the actual 
conditions (proximity to phase transition points and the typ- 
ical size of emerging spatial patterns) the linear system size 
was varied from L = 200 to 1600 and the relaxation time was 
varied from 10 4 to 10 7 MCS to ensure proper accuracy. 

In Fig. Q] we plot the critical value of r, defined as the min- 
imally required R = r/G (note that the normalization with 
G enables comparisons with results on graphs other than the 
square lattice 13011 ') where pc first becomes > 0, in depen- 
dence on A. It can be observed that the stronger the Matthew 
effect, the lower the multiplication factor needed for the sus- 
tenance of cooperation. This holds irrespective of K. While 



FIG. 2: (Color online) Critical R at the D — >• C + D phase transi- 
tion in dependence on A for three different interaction graphs. In- 
sets show the graphs schematically in the color (symbol type) cor- 
responding to the results in the main panel, with the group size G 
indicated for vertices encircled red (gray). As in Fig.fT]the promo- 
tion of cooperation is clearly inferable, irrespective of the properties 
of the interaction graph. Here K = 0.1 was used in all three cases. 
The dashed orange line shows the result as predicted by the mean- 
field approximation (see main text for details). 



A = 0.001 essentially returns the classical version of the 
game, with the same R as reported in [30], a sharp descent 
towards as low as R = 0.17 follows for larger A. The ef- 
fect saturates and may even revert slightly when A becomes 
comparable to G, which is due to the fact that M x must not ex- 
ceed G since all the accumulated public goods within a group 
are then already assigned to player x. To test the generality 
of the observed promotion of cooperation further, we plot in 
Fig.|2]the critical R for three different interaction graphs (see 
insets) with very distinct properties that are know to vitally 
affect the evolution of cooperation [31]. Namely, a regular 
graph with zero clustering coefficient (square lattice), a regu- 
lar graph with a high clustering coefficient (triangular lattice), 
and a random regular graph having no local structure. Irre- 
spective of these details the Matthew effect promotes cooper- 
ation equally well, as can be concluded from the descending 
i?(A) phase transition lines. 

An inspection of the stationary distribution of M x values 
for intermediate A reveals a double peaked Gaussian with 
maxima approximately at — G/2 and G/2 respectively. This 
is easily traced back to the nature of the Matthew effect 
in that it segregates the population into successful (peaked 
around G/2) and unsuccessful (peaked around —G/2) play- 
ers. What is more intriguing is the observed promotion of 
cooperation, which was in the past typically associated with 
strongly heterogeneous (e.g. power law or exponential) dis- 
tributions, whereby it was argued that the high-impact (in 
our case equivalent to high M x ) players stabilize cooperation 
while defectors are doomed by means of a negative feedback 
effect [ 31|] . This well-known explanation applies only par- 
tially in our case. In fact, here cooperators benefit also from 
the dynamical reshaping of the evolutionary landscape (de- 
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FIG. 3: (Color online) Strategy regions in the A — R parameter space 
for the square lattice at K — 0.1. While R needed for cooperation 
to survive decreases steadily with increasing A (solid red line), the 
opposite holds for the extinction of defectors. The transition to the 
pure C region (dotted green line) is preceded by the emergence of a 
peculiar C +D region (dashed blue line), where a minute fraction of 
defectors is able to survive in the presence of dominating coopera- 
tors. The inset shows two characteristic time courses for just above 
[dashed green line (top arrow at R — 1.5)] and below [solid blue 
line (bottom arrow at 7? = 1.47)] the C +D — > C transition line 
for A = 0.001. Within the C +D region defectors utilize their one- 
in-a-million survival chance to evade extinction. Although in the 
inset only the first 4096 MCS are depicted, we have verified that at 
R = 1.47 pr> > for up to 10 7 MCS at L = 1600 system size. The 
dashed black line in the inset is proportional to MCS . 



fined by M x ), which renders the defective strategy maladap- 
tive and thus gives way to a mixed C + D phase. This can 
be confirmed analytically by means of a simple well-mixed 
approximation of our model. Namely, by plugging into the 
well-mixed ansatz i32ll p'c = pcO- ~ Pc){Pc — Pd) the pay- 
offs as defined above, treating Nc/G as pc and introducing 
Gaussian distributed values £ with standard deviation a for 
the difference in M (as motivated by the observed distribu- 
tion of M x values), we obtain p'c = PcX}- ~ PcK r Pc£ — !)■ 
Following a first-order small-noise expansion 13311 we finally 
have p c = pc(l - Pc)(v 2 rpc{pc ~ ZPc/" 2 ) ~ !)■ From 
p'c = and the sign of the second derivative we find pc = 
and pc = 1 as unstable steady states. The new stable steady 
state of pc comes from <J 2 rpc{pc — 3p c /2) — 1 (the poly- 
nomial has three roots, two of which complex conjugates), 
which however is to cumbersome to be expressed here explic- 



itly. What is relevant is that one obtains r oc cr~ a5 for the 
D — s- C + D phase transition, which can be verified easily by 
straightforward numerical integration of p'c- Since we lost G 
in the well-mixed approximation r may be rescaled to account 
for R as obtained from Monte Carlo simulations and a takes 
on the role of A. Preserving the slope of the dependence this 
yields the dashed orange line depicted in Fig. [2] While it is 
by no means implied that this approximation gives a good fit 
to the Monte Carlo simulations, it nevertheless confirms the 
effect by means of a simple analytically treatable model. 
The promotion of cooperation, however, is only one ob- 
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FIG. 4: (Color online) Uppermost panel shows time courses of 
M x for the "super-persistent" defector (solid black), its four near- 
est neighbors (dashed red), its four next-nearest neighbors (dash- 
dotted green) and the four next-next-nearest neighbors (dotted blue). 
The three middle panels depict gray-scale coded (white maximal, 
black minimal) values of M x around the "super-persistent" defec- 
tor at times corresponding to the vertical dashed lines in the upper- 
most panel. Colored stars serve the identification of the neighbors 
corresponding to the colors (shades of gray) used in the uppermost 
panel. The bottommost row depicts possible configurations of de- 
fectors (gray) around the "super-persistent" defector (black). Re- 
sults were obtained on a square lattice of linear size L — 1600 at 
K = 0.1, R = 1.47 and A = 0.001. 



servation following from the introduction of the Matthew ef- 
fect. There is another, namely the preservation of defectors 
to be demonstrated below, which is rooted in extremely rare 
small-region effects originating from the spatiality of the con- 
sidered model. Figure [3] depicts the strategy regions in the 
A — R parameter space for the square lattice, from where the 
peculiar behavior can be inferred. In particular, as the crit- 
ical R inducing D — » C + D falls, the critical R required 
for C + D — > C raises sharply with increasing A. Extensive 
Monte Carlo simulations using L = 1600 reveal a narrow in- 
termediate region denoted as C +D , where the most minute 
fraction of defectors can prevail indefinitely in the sea of co- 
operators. The inset features two time courses just above and 
below the C +D — > C transition line (see arrows in the main 
panel). What at first appears to be algebraically slow relax- 
ation (marked by the dashed line in the inset) is facilitated by 
the coevolutionary impact of the Matthew effect to become ei- 
ther an absorbing C region (green) or the C +D region (blue), 
where approximately one defector is able to survive amongst 
10 6 cooperators (note that the temporal courses are averages 
over 10 independent realizations). We also note that we use 
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region rather than phase for the C + D — > C +D — > C tran- 
sitions since we were unable to completely resolve the large 
system size limit. 

Zooming in on the neighborhoods around these "super- 
persistent" defectors gives vital clues with regards to their sur- 
vival. Figure [4] shows the time evolution of M x for one such 
defector and its neighbors. Evidently, it is possible, although 
highly unlikely, that isolated defectors become fully protected 
by their neighbors by means of the spontaneous emergence 
of successive microscopic hierarchies (expressed in terms of 
M x ) in space. The three middle panels show that while the 
central player (defector) enjoys the benefits of the highest pos- 
sible M x , its four nearest neighbors have the lowest possible 
M x associated to them. Further crucial is then the fact that 
all the neighbors of the four nearest neighbors, like the central 
defector, also all eventually acquire the highest possible M x . 
This makes it impossible for the four potential donors of the 
new strategy to the central defector to overtake it. Even though 
they may occasionally succeed in adopting the cooperative 
strategy, it is impossible to enforce the latter onto the central 
defector because they are at an inherent disadvantage due to 
the unfavorable outcome of the Matthew effect in their imme- 
diate neighborhood. Consequently, the only possible configu- 
rations (showing defectors only) around the "super-persistent" 
defector are those depicted in the bottom row of Fig. 0] with 
zero chances of full cooperator dominance. It is to be em- 
phasized that although one defector amongst a million co- 
operators seems negligible, from the evolutionary viewpoint 
it nevertheless preserves the seed for defection to take on a 
more dominant role should the conditions ever become more 
favorable. From the view point of statistical mechanics, it is 
fascinating to learn that coevolutionary rules may, due to the 
interplay of spatial structure and the dynamics of the coevolv- 
ing quantity, induce spatially extremely localized microscopic 
patterns that prevent phase transitions into absorbing states 
despite the fact that the strategy destined to die out is com- 
pletely maladaptive (note that in the absence of the Matthew 
effect cooperators reach full dominance on the square lattice 
at i? w 1.1 10]). 

In sum, we have demonstrated that a simple strategy- 
independent coevolutionary rule mimicking the Matthew ef- 
fect promotes cooperation in group interactions amongst unre- 
lated individuals. This can be partially attributed to the emer- 
gent Gaussian distribution of factors determining the distri- 
bution of public goods, but also to the inability of defectors 
(or their lesser ability if compared to cooperators) to adapt 
to continuously changing evolutionary landscapes, as demon- 
strated by an analytically treatable well-mixed approximation. 
However, the Matthew effect may also give rise to spatially 
highly localized microscopic patterns that protect maladaptive 
strategies despite of their obvious evolutionary disadvantage, 
thus preventing absorbing phases and preserving an option of 
a comeback of seemingly exterminated traits. We hope that 
this work will inspire further research in the yet unexplored 



directions concerning the impact of coevolutionary rules on 
the outcome of games governed by group interactions. 
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